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ABSTRACT 

Binaries that contain a hot subdwarf (sdB) star and a main sequence companion may have interacted 
in the past. This binary population has historically helped determine our understanding of binary 
stellar evolution. We have computed a grid of binary population synthesis models using different 
assumptions about the minimum core mass for helium ignition, the envelope binding energy, the 
common envelope ejection efficiency, the amount of mass and angular momentum lost during stable 
mass transfer, and the criteria for stable mass transfer on the red giant branch and in the Hertzsprung 
gap. These parameters separately and together can significantly change the entire predicted population 
of sdBs. Nonetheless, several different parameter sets can reproduce the observed subpopulation of 
sdB + white dwarf and sdB + M dwarf binaries, which has been used to constrain these parameters 
in previous studies. The period distribution of sdB + early F dwarf binaries offers a better test of 
different mass transfer scenarios for stars that fill their Roche lobes on the red giant branch. 
Subject headings: binaries: close — subdwarfs — stars: horizontal branch 



1. INTRODUCTION 

Binary population synthesis (BPS) models parameter- 
ize several physical processes that can drastically alter 
man y properties of th e stellar populations that they pre- 
dict ( jEggleton 2006). While some of these parameter- 
izations are used to reduce the computational burden 
of evolving > 10 6 binaries, many stem from our incom- 
plete understanding of single and binary star evolution. 
These phenomenological parameterizations ar e tuned to 
reproduce w ell-known binary populations (e.g.|Nelemans| 



son 



1990[ lYungelson fc Tutukov||2005[ |Nelemans fOTu| 



tion with a theoretical approach (e.g. |Tutukov k, Yungel 

)05~ 



Han et al.l (|2002p computed detailed binary evolution 



models for five sdB formation channels: the first and 
second common-envelope (CE) ejection channel, the first 
and second stable Roche lobe overflow (RLOF) channel, 
and the double helium white dwarf (WD) merger chan- 
nel. The results of these cal culations w e re the n applied 

hereafter 

H03) 



in a BPS study presented in Han et al 

able 



J2003 
to reproduce certain 



et al. 2000fo. Binaries containing a subdwarf B (sdB) 



star constitute one often-used population. These objects, 
also called extreme horizontal branch (EHB) stars, are 
thought to be core helium burning stars with very thin 
hydrogen envelopes, M env < 0.02 Mq, (for a recent re- 
view, see |Heber||2009| . Such objects can be formed in 
interacting binaries, by stripping the hydrogen e nvelope 



The BPS model was 
properties of the observed population, including the bi 
nary fraction and the mass es and periods of the sdB 



+ WD binaries reported i n [Maxted et al. (|2001|) and 



|Morales-Rueda et al.| (]2003|) 

aspects of this mode 



)espite the success of many 



from a star as it ascends the red giant branch (Mengel 
Observations support a binary origin for a 



et al. 



1976) 



significant fraction of sdBs. Not only are mos t known 
sdBs in bin aries with tight periods P < 10 days ( |Maxted| 
et al ||2001| , but their semimajor axis distribution is m- 
consistent with the primord i al distribution prop osed by 
Duquennoy fc Mayor| (|1991 [) ( |Heber et aL"1|2002j) . Being 
ubiquitous, bright, and produced by interacting binaries 
hot subdwarfs are an ideal population to test binary- 



evolution scena rios (|Green et al. 
Nelemans|[2010 ). 



2001 Han et al. 2003 



Several previous studies have investigated sdB forma- 



some properties of the modeled 
population are not consiste nt with a recent sdB radia l 
velocity study discussed in Copperwheat et al. (2011). 
These authors found that the binary fraction among sdBs 
is 50%, but it is the absence of short period (P < 10 days) 
sdB + G or K type dwarf binaries from their sample that 
is most incompatible with the predictions of H03. For the 
sdBs t hat showed evidence of a G or K type companion 
in the Copperwheat et al. (20111 sample, no radial ve- 
locity variations were detected on short (few day) time 
scales, which implies that these systems have long peri- 
ods. They concluded that the presence of a G or K type 
companion is indicative of a long period binary, which is 
co nsistent with the s dB radial velocity survey carried out 
by Green et al. (2001 ) who inferred that P ~ 1 — 3 years 
for such systems. Only one short period sdB + G or K 



Electronic address: dclausen@astro.psu.edu, wade@astro.psu.edu, 
ravi@gravity.psu.edu, oshaughn@gravity.phys.uwm.edu 



dwarf system has been reported, Moni Bidin &: Piotto 
(20101 identified such a binary in the globular cluster 
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NGC6752 with P < 5 days. This is in stark contrast to 
the large population of sdB + G or K type systems with 
P < 10 days predicted by H03, see, e.g., their Figure 15. 

To explore the discrepancy between the models and ob- 
servations and to test how the assumptions made in bi- 
nary stellar evolution models impact the sdB population, 
we have carried out our own BPS study of sdB forma- 
tion. In lscction 21 we describe our BPS ca lculation. The 
results of our calculation are discussed in Isection 3l We 
show the initial binary distributions that lead to sdBs 
with hydrogen-burning companions; we show the joint 
distribution of orbital periods and companion effective 
temperatures at the present epoch; and we discuss which 
combinations of BPS parameters are compatible with ob- 
servations of sdB + WD and sdB + M dwarf populations. 
We describe how sdB + early F binaries can better con- 
strain the models in Isection 4l In Isection 5l we com- 
pare our models with previous results. Throughout, we 
distinguish between true helium-burning sdB stars and 
post-RGB stars with similar effective temperatures and 
surface gravities. 

2. BINARY POPULATION SYNTHESIS 

In this work, we focus on binaries in which the pri- 
mary, the initially more massive star, has evolved off of 
the main sequence (MS), expanded to fill its Roche lobe, 
and lost its envelope to become an sdB. The envelope 
can be lost either through stable RLOF mass transfer 
to the secondary or through CE ejection. In these sys- 
tems, the secondaries are typically MS stars, but some 
will have evolved as far as the red giant branch (RGB). 
For simplicity we will often refer to all binaries in which 
the sdB has a hydrogen burning companion (both core- 
burning MS stars and shell-burning RGB stars) as sdB 
+ MS binaries. Although we were interested in the pop- 
ulation of sdBs with hydrogen burning companions, we 
also tracked sdB + WD binaries, in which the secondary 
has evolved to become the sdB, because nearly all of the 
observational constraints come from this subpopulation. 

2.1. Binary Evolution Code 

We used the binary evolu tion code BINARY STAR 
EVOLUTION (BSE) described in |Hurley et~aL] ( |2002h . Be- 
low we describe several aspects of BSE that are r elevant 
to our study , but we refer the interested reader to [Hurley] 



et al. (20021 for a complete description of the code. In 
our models, we varied several parameters that govern a 
system's evolution to explore how each of these compo- 
nents impact ed the re sultant population of sdB + MS 
binaries (see Table 1). First, we varied the minimum 
core mass at which helium ignites. By default in BSE, 
stars that undergo a helium flash (effectively stars with 
Mzams < 1-995 M at solar metallicity) do not ignite 
helium if their envelopes are lost before they reach the 
tip of the RGB (TRGB). However, [P7Cruz et al.| ( [19961 ) 
showed that it is possible for these stars to ignite helium 
if the envelope is lost within ~ 0.4 bolometric magni- 
tudes of the TRGB. Stellar evolution models described 



Han et al. (2002 1 suggest that if the mass loss is rapid, 



as is the case for CE ejection, stars can ignite helium if 
the core has reached 95% of mass it would have at the 
TRGB (m c (TRGB)) before the envelope is removed. We 
modified BSE to allow stars to ignite helium if the core 



has reached a mass 



/He TOc (TRGB) 



(1) 



b efore the star's en velope is lost , accordi ng to the result 
of Han et al. ( 2002 1 . As seen in Table 1 we constructed 
synthetic populations with /ho 0.95 that allowed he- 
lium ignition in stars that lose their envelopes shy of the 
TRGB, through either CE ejection or stable RLOF, and 
models with /h g = 1 that required the stars to reach the 
TRGB before helium ignition. 

Second, when the donor star fills its Roche lobe, the 
mass ratio q = M^/M^ determines whether or not mass 
transfer is dynamically stable. Here Md is the mass of the 
donor star and M a is the mass of the accretor. Stability 
of mass transfer depends on the structure of the donor, 
the mass ratio of the binary, and other factors like the 
amo unt of mass and angular momentum lost to infin- 
ity ([Hjellming fc Wcbbink 1987} |Soberman et al.||1997 



Wcbbink 2006; |Ge et al.||2010p When q is smaller than 
a critical value g cr it the mass transfer is stable, otherwise 
the system will undergo CE evolution. The functional 
form of the critical mass ratio versus all stellar and bi- 
nary parameters is poorly known. To gauge how differ- 
ent assumptions about g cr it influenced the population of 
sdB + MS binaries, we considered many different values 
for this parameter. For systems undergoing RLOF while 
the donor is in th e Hertzsprung gap (HG), we set q cr it 



to either 3.2 or 4 ( |Han et al.||2003| |Hurley et al. 2002 
respectively). If the system undergoes RLOF while the 
primary is on the RGB, BSE's default critical mass ratio 
is given by 



q alt = 0.362 



1 



3(1-^) 



(2) 



where m c d is the donor's core mass and Md is the total 
mass of t he donor sta r. We have used the value of g cr it 
given by |Equation 2| in some of our BPS models. Un- 
der this assumption, q CT i t < 1 in most cases and mass 
transfer is unstable. |Equation 2 was derived from mod- 
els of conservative mass transfer between two condensed 



polyt r opes b y Hjellming & Webbink ( 1987 ) , however Han 
et al. ( 2002 ) described a set of detailed binary stellar evo- 
lution calculations for completely non-conservative mass 
transfer in which the mass lost from the system carries 
away the specific angular momentum of the accretor. Un- 
der these conditions, mass transfer was stable for q < 1.2. 
H03 argued that enhanced wind mass loss would have 
the effect of increasing q cr it- Following this suggestion, 
we also computed models with q cr - lt — 1.5 on the RGB. 

Third, in addition to changing the criteria that deter- 
mine whether or not mass transfer will be dynamically 
stable, we varied parameters that govern each mass loss 
scenario. During stable RLOF mass transfer, some of 
the mass lost by the donor might not be accreted by 
the companion, but instead will be lost from the system. 
We defined /rlof as the fraction of mass lost by the 
primary that is accreted by the secondary during stable 
RLOF mass transfer, e.g., /rlof = 1 corresponds to con- 
servative mass transfer. By default in BSE, the amount 
of material accreted by a MS or HG secondary is limited 
by its thermal timescale and /rlof is computed with 



/rlof = min 1.0, 10 



M a 



M d T K H 8 



(3) 
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where M a is the mass of the accretor, M<j is the mass 
loss rate of the donor, and tkhe is the Kelvi n-Hclmholtz 
timescale of the accretor. Setting /rlof with |Equation 3] 
allows the accretor to accept all of the material lost by 
the donor for sufficiently low M^, but caps the accretion 
rate at a value that prevents the accretor's radius from 
expanding by more than a factor of ~ 1 .5 for larger Md 
( |Tout et al.|ll997| JPols & Marinus||1994[ ). We also com- 
puted models with /rlof = 0.5, again, following H03. 

Fourth, the mass lost by the primary that is not ac- 
creted by the secondary leaves the system and it must 
also carry away angular momentum. This material can 
take angular momentum from the orbit, the donor, or the 
accr etor. In the first case, denoted by 7 = 1 in the text 
and |Tablc l[ the material takes away a specific angular 
momentum proportional to the orbital angular momen- 
tum of the system, a 2 £!orb where fiob is the orbital 
angular frequency and a is the semimajor axis. In the 
second case, denoted by 7 = — 1, the material takes away 
a specific angular momentum proportional to that of the 
donor, a^fJorb- Cases in which the material carries away 
a specific angular momentum proportional to that of the 
accretor, a^orb, are denoted by 7 = —2. Note that q cr it 
and the parameters that describe non-conservative mass 
transfer, /rlof arL d 7, are coupled for fixed assumptions 
about stellar structure, but in our parameter study we 
treat them as independent inputs. 

Finally, we also considered several different values for 
the parameters that govern CE evolution. During a CE 
phase, the mass donor's envelope drastically expands and 
engulfs the accretor. While this process is not well un- 
derstood, it is believed that the core of the donor and 
the entire accretor spiral towards one another within the 
CE as orbital energy goes into heating and unbinding 
the CE. This process drastically reduces the binary's or- 
bital period, and often results in a merger. In BSE this 
process is controlled by two parameters. The struc- 
ture parameter A determines the envelope binding energy 
and cuce determines how efficiently the orbital energy is 
transferred to the envelope. We considered two cases 
for the envelope binding energy. One was the BSE de- 
fault, which is to determine the value of A from the 
star's structure and evolutionary state using fitting for- 
mulae derived from single star evolution models. With 
the value of A the binding energy (BE) is then calculated 
as Sbind = GMi(Mi — m c i)/(Ai?i), where Ri is the ra- 
dius of the mass donor. Alternatively, we used the ana- 
lytic expressions for th e envelope binding energy given in 
Loveridge et al. (2011 hereafter LVK11). These authors 
computed the envelope binding energy from a grid of de- 
tailed stellar evolution models and fit it as a function of 
metallically, mass, radius, and evolutionary phase. The 
fits we used included terms for the internal energy of the 
gas, the radiation energy, and the gravitational energy, 
but not the recombination energy. As the fraction of the 
orbital energy available to eject the common envelope, on 
physical grounds «ce < 1- Treating the CE parameter 
qce as a purely phenomenological factor, we considered 
«ce = 0.75, 1.5, and 3. 

2.2. Initial and Present Day Population 

The initial population of main sequence binary sys- 
tems, some of which eventually have an sdB companion, 



was chosen as follows. To ensure that the distribution of 
sdB binaries is well sampled, we selected 3.5 x 10 6 ini- 
tial binary systems, with the mass of both the primary 
(-M10, more massive) and secondary (M20, less massive) 
in the range 0.1 Mq — 10 Mq. The primary's mass dis- 
tri bution was taken from the initial mass f unction (IMF) 
of |Kroupa fc Weidner (2003), as given in Hurley et al. 
(2002) 



£(m) 



0, 




to < mo 


bim~ 


1.3 


to < to < 0.5 




2.2 


0.5 < to < 1.0 


bim~ 


2.7 


1.0 < TO < OO 



(4) 



where t;(m)dm is the probability that a star has a mass 
between to and m + dm and we have used too = 0.1 Mq. 
The initial mass ratio (l/<Zo) distribution was taken to 
be uniform between and 1, so the secondary's mass 
was correlated with the primary by M2 — Mi /go- The 
initial mass ratio distribution is not well constrained by 
observations, but has a significant impact on the results. 
For example, H03 considers both a constant mass ratio 
distribution case and a case in which each member of 
the binary is chosen form the IMF independently. In the 
uncorrelated mass distribution case, the number of sys- 
tems formed through CE ejection increased drastically. 
However, an uncorrelated m ass ratio distribution is dis - 
favored by observations (e.g. Duquennoy fc Mayor|1991 ). 
We did not systematically explore the mass ratio distri- 
bution, and chose to use the constant mass ratio distribu- 
tion in all of our runs to enable comparisons with H03's 
best fit model, which also used a uniform distribution 
of 1/qo- We note, however, it is possible to rescale any 
simulation to an arbitrary mass ratio distribution. The 
distribution of initial orbital separation (a) was t aken to 
be uniform in ln(a) between 3 Rq and 10 4 Rq (Hurley 
et al. |2002 , and references therein) . All orbits were as 
sumed to be circular and all stars were assumed to have 
solar metallicity. 

The same initial population of main sequence binaries 
was given as an input to the BSE code for each run. The 
code uses an adaptive integrator to sample the evolution- 
ary tracks of each binary. In order to obtain the present 
day distribution of sdB binary systems produced during 
the evolution of these initial main sequence binaries, we 
resampled the output to a fixed sampling rate, linearly 
interpolating the characteristics of each individual star 
(gravity, effective temperature, luminosity etc.) so that 
we have a datapoint for every 10 Myr, starting from its 
birth, until either the maximum evolution time (15 Gyr 
in our case) or until one or both of the stars had no rem- 
nants. We then selected systems in which one or both 
the stars are sdBs. 

Since we cannot observe whether a star is burning 
helium in its core, we chose to identify sdBs by their 
positions in (T c ff,logg) space. For this work we de- 
fined the "hot subdwarf box" to be the region where 
20000 < T off < 45000 K and 5.0 < logg < 6.6 (cgs). 
Once a star enters the "hot subdwarf box" we record the 
resampled evolutionary history, while the system remains 
in an sdB phase. The collective evolutionary record of 
all such systems that are in the sdB phase then mimics 
a sample of present day sdB binaries that are in various 
stages of evolution. Using this method, the lifetimes of 
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Fig. 1. — Cartoon of the regions of log P — log T c g plane occupied 
by sdB binaries produced by different evolutionary channels. Here 
log T c ff refers to the temperature of the MS companion. Note that 
these regions are the composite of all of our runs. With a particular 
set of parameters, only a subset of these regions might appear and 
only a subregion of each region may be populated. 

systems in an sdB phase are appropriately reflected in 
the present day population. If an sdB binary has a short 
lifetime, it remains in the sdb "box" for fewer resampled 
timesteps and hence such a system contributes less to 
the total number of current systems. However, we do 
not specify a star formation rate and defer discussion of 
the birthrates or space densities of the binaries produced 
in our BPS models to future work. 

Furthermore, we only consider sdBs with companions. 
There are many proposed formation channels for single 
sdBs, including the merger of two He white dw arfs (WDs) 



(Webbink 1984 Iben k Tutukov||1984 1986), enhanced 



RGB mass loss ( jD'Cruz eT al. 1996), ejection of the H 



envelope by a sub-stellar companion (Soker 1998), cen 
trifugally enhanced mass loss triggered b y common enve 
lope (CE) mergers (jPolitano et al.|2008|), and the merge r 



of an M dwarf with a He WD (Clausen k Wade 
However, the algorithms in BSE either do not inc 



2011] 
ude or 



do not self-consistcntly model these mergers and the evo- 
lution of the remnant. We therefore draw no conclusions 
about the fraction of sdBs that are in binaries. 

We note that the mass transfer and stellar evolution al- 
gorithms used in BSE do not produce true sdBs, instead 
they produce naked helium stars that have had the entire 
hydrogen envelope removed. These naked helium stars 
all fall along the He main sequence, which runs through 
the high log fl, high T e g c orner of our sdB box. Using 
the models of Caloi ( 1972[ ), we have confirmed that these 
naked helium stars would still be in the "hot subdwarf 
box" if they had hydrogen envelopes of up to 0.02 Mq. In 
the most extreme case, adding a 0.02 Mq envelope to the 
least massive star in our models (M = 0.32 Mq), logg 
would decrease by 1.2 dex, logT ff would decrease by 0.2 
dex, and the star would still be in the "hot subdwarf 
box." Although BSE does not produce sdBs with hydro- 
gen envelopes, we have only selected systems that would 
still be in the proper log g and T e s range if they had 
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Fig. 2. — Cartoon of regions of the log Pq — Mio plane that pro- 
duce helium-burning sdB + MS binaries through four evolutionary 
channels and the region that produces post-RGB stars that sit in 
the sdB "box". 

a thin hydrogen envelope. Because of this limitation of 
BSE, we do not discuss the T e g or log g distributions of the 
sdB stars themselves. Nor without further assumptions 
can we discuss the radiative properties (magnitudes, col- 
ors) of the sdBs or the binaries in which they are found. 
The sdB masses (with small errors from neglecting the 
H envelope mass) remain valid for discussion, however. 

3. RESULTS 

The sdB binary popul ations pr oduced in our BPS mod- 



els are summarized in Table 1 Column (1) gives the 



identification number for the run and the next seven 
columns list the parameters and parameterizations used 
in the run. Column (9) gives the number of initial bina- 
ries that produce an object that lies in the sdB box and 
column (10) lists the number of resampled, present day 
binary systems containing such an object. Column (16) 
gives the percentage of the resampled systems that con- 
sist of a core helium burning sdB with a hydrogen burn- 
ing companion and the preceding five columns give the 
percentage of systems formed by each formation channel. 
Column (17) lists the percentage of core helium burning 
sdB binaries with WD companions. Column (18) lists 
the percentage of present day systems that do not con- 
tain a true helium-burning sdB but instead an post-RGB 
star that is presently in the sdB box on its way to be- 
coming a W D. The model suitability parameter p (see 
section 3.1) is listed in column (19). 

Our BPS models show that the population of sdB + 
hydrogen burning companion binaries is highly sensitive 
to the value of several parameters necessary to model bi- 
nary evolution. Varying the values acE, 7> <?crit, /rlof, 
/hg, and the envelope binding energy, within reasonable 
ranges for each, resulted in substantial changes to the 
distributions of both the orbital periods and companion 
masses for sdB + MS binaries. Furthermore, our models 
suggest that these parameters are poorly constrained by 
the current sample of sdB + WD and sdB + M dwarf 
binaries with measured orbital periods because many of 
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our models produced a population of such systems that 
was consistent with the observed sample. Assumptions 
about the IMF, the distribution of initial orbital peri- 
ods and eccentricities, and the distribution of the initial 
mass ratios also affect the the populations and our con- 
clusions, but were not varied among the different models 
presented. 

3.1. General Description 

Before we discuss and compare the populations pro- 
duced by individual runs in detail, we present a gen- 
eral description of the results and outline how we will 
present, examine, and evaluate each run. To display the 
observable properties of these predicted populations, we 
will plot each binary's orbital period against the effective 
temperature ( T c g) of the MS companion in logP— logT c ff 
diagrams like |Figurc lj The black lines outline the re- 
gion occupied by the sdB + MS binaries predicted by the 
best-fitting BPS model chosen in H03 (see their Figure 
15b). Also shown are sketches of the regions occupied by 
sdB + MS binaries formed through different formation 
channels if we combine all of our BPS models. The re- 
gions shown are schematic; the exact boundaries of these 
regions and whether or not they are populated for a par- 
ticular run depend on the parameters used in that run. 

Systems that form as the result of stable mass transfer 
beginning while the primary is in the HG (HG RLOF) 
tend to have the most massive companions because up to 
3.5 Mq of material is transferred to the companion. How- 
ever, the period distribution of these systems is sensitive 
to assumptions about the amount (/rlof) and specific 
angular momentum (7) of material lost to during mass 
transfer. 

If we set cfo-it = 1.5 on the RGB, then sdB + MS bina- 
ries can form through stable mass transfer that begins on 
the RGB (RGB RLOF). In these binaries, the MS com- 
panions span a large range in temperature. This is due, 
in part, to evolution of the companion towards the RGB 
while the primary is an sdB. The periods of these systems 
are typically longer than 5 days. Some binaries survive 
CE evolution that begins on the RGB (RGB CE) to form 
sdB + MS binaries. In the vast majority of these sys- 
tems, the companion is of type G or later. Note that the 
H03 models presented predict systems with companions 
to a lower T c g than our models do. This is a consequence 
of our choosing too = 0.1. If we instead take too = 0.08, 
then the log T c g range extends to 3.33. Since we do not 
draw M 2 from the IMF, but determine it as M 2 — qoMi, 
there are few systems with companion masses this low. 

We selected sdBs from the BPS runs based on their 
surface gravity and effective temperature, much as an 
observer would. By this process, we also selected many 
objects that were not core helium burning sdBs, but the 
cores of post-RGB stars cutting across the "hot subdwarf 
box." These objects are less massive than helium burn- 
ing sdBs and have an average mass of 0.2 Mq. Individ- 
ually, these systems spend little time in the sdB "box," 
however they are quite numerous. Interestingly, many 
of these post-RGB binaries have B-F type companions. 
These systems have similar periods and WD masses to 
four early type dw arf + He WD systems recently discov- 



them in the plots below. However, for clarity we will 
refer to them as post-RGB stars and any discussion of 
sdBs below refers only to core-helium burning sdBs. 
Figure 2 shows the initial primary masses Mm and 



orbital periods Po of the binaries. The systems that form 
sdBs through the channels labeled in |Figurc~T1 come from 
distinct regions in the log P — M w plane. Again, these 
regions are schematic and represent the combination of 
all of our BPS runs with the boundaries changing slightly 
depending on the parameters used. The lower limit on 
M10 for sdB progenitors is determined by the 15 Gyr 
evolution time used in our BPS models. Had we used a 
10 Gyr evolution time, the limit on M w would shift from 
~ 0.8 Mq to ~ 1 Mq. Figure 2 shows this would reduce 



Breton et al. 2011D. These objects meet the observable- 



the number of sdB binaries produced by the RGB CE 
and the RGB RLOF channels. 

Systems that evolve via the HG RLOF channel have 
primaries with M\q > 2 Mq and orbital periods Pq < 10 
days. Binaries with initial periods below the lower bound 
shown in |Figurc 2| begin RLOF while the primary in on 
the MS, which results in either the primary losing too 
much mass to ignite helium in its core or a merger of the 
two stars. The post-RGB systems undergo stable RLOF 
mass transfer that begins while the primary is in the 
HG, but the envelope is stripped before the core reaches 
the mass required for helium ignition. Accordingly, these 
systems occupy the low-mass, short-period corner of the 
diagram. For binaries with longer Po, the primary does 
not come into contact with its Roche lobe until it begins 
to ascend the RGB. If we allow g cr it = 1-5 on the RGB, 
some of these systems can undergo stable mass transfer 
and form sdBs through the RGB RLOF channel. At 
longer periods still, the star will ignite helium and detach 
from its Roche lobe before the envelope is completely 
stripped. This forms the long period boundary for RGB 
RLOF systems. 

If we use the BSE q c[ - lt formulation, then systems that 
begin mass transfer while the primary is on the RGB 
will enter a CE phase. We have split this channel into 
two regions, RGB CE I and RGB CE II. In the for- 
mer, M w < 1.995 Mq and these stars only produce 
sdBs if they can ignite helium below the TRGB. The pri- 
maries in RGB CE II ignite helium non-degenerately on 
the giant branch. The lower boundary of the combined 
RGB CE regions is approximately given by Pi (TRGB) = 
Pl(Po)9o = 1), where Pi is the radius of the primary 
and Pl is the Roche lobe radius. Binaries with smaller 
go occupy the region above the curve because we have 
included tidal evolution in all of our models. Without 
tidal evolution, the Pi (TRGB) = Pl(Po,(?o = 1) curve 
would be the upper boundary of this region. Tidal evo- 
lution, which in BSE drives the binary to a shorter orbital 
period through spin-orbit angular momentum exchange, 
is stronger in systems with large qo ( Hut|198l ). Without 
tidal evolution a binary with M\q ^ 1.995 and qo < 1 
would need a shorter Pq for the primary to fill its Roche 
lobe near the TRGB. We note that the upper boundary 
of the RGB CE I region in our models is flat because we 
require M 2 > 0.1 M . 

To evaluate each set of parameters we compared the 
modeled populations with the observed properties o f sdB 



ered with Kepler (IRowe et al. 2010 Carter et al. 2011 binaries presented in |Morales-Rueda et al.| p003) and 



Copperwheat et al. ( 201ip . To replicate selection biases 
based criteria we have used to select sdBs, so we include in the observed sample, we only compare the sdB + WD 
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Fig. 3. — logP vs logT c ff of the MS companion for 6 BSE runs. In each panel, 7=1 and c*ce = 0.75. Blue crosses and red asterisks 
represent systems that underwent stable mass transfer beginning while the primary was in the HG or RGB, respectively. Systems that 
underwent CE evolution that began while the primary is on the RGB are plotted with green diamonds. Binaries that contain a post-RGB 
object instead of a core helium burning sdB are plotted as purple circles. The black lines give a rough outline of the region occupied by 
systems from H03. Panel (a) shows our fiducial run using all BSE default values, Run 1. Panel (b) shows Run 2, identical to Run 1 but 
setting /ho = 0.95. Panel (c) shows Run 3, identical to run 2 but using the LVK11 envelope binding energies. Panel (d) shows Run 4, 
identical to Run 3, but setting g cr it = 1.5 on the RGB. Panel (e) shows Run 5, identical to Run 4 but setting /rlof = 0.5. Panel (f) 
shows Run 6, the same as Run 5 but with the BSE envelope binding energies. The temperature regimes for spectral types B-M are shown 
in panel (a). 
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Fig. 4. — logPo vs Mio (initial orbital period and primary mass, respectively) for the BSE runs shown in |Figure 3| Crosses and asterisks 
represent systems that will undergo stable mass transfer beginning while the primary is in the HG or RGB, respectively. Systems that will 
undergo CE evolution that begins while the primary is on the RGB are plotted with diamonds. The circles are systems that produce a 
post-RGB object instead of a core helium burning sdB. The color of each point shows the system's initial mass ratio qo (see color bar at 
top of figure), (a) Run 1, (b) Run 2, (c) Run 3, (d) Run 4, (e) Run 5, and (f) Run 6. 
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Fig. 5. — Simulated minimum companion mass (see |section 3.2[l 
diamonds) and sdB + M dwarf (green crosses) binaries are plotted. 

binary "obse rved " at different orientation a ngles. The red symbols are an observed sample of sdB binaries compiled from Moralcs-Rucda 
|et al.| |2003| l and Coppcrwh eat et al.|l|2011| . The squares are sdBs with WD companions, the circles are sdBs with M dwarf companions, 
and the triangles are systems with unknown companion type, (a) Run 1, (b) Run 2, (c) Run 3, (d) Run 4, (e) Run 5, and (f) Run 6. 
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Fig. 6. — Same as |Figure 3| with, additionally, light blue triangles representing systems that underwent CE evolution beginning while 
the primary was in the HG tor (a) Run 17 showing the effect of setting ciqe = 1-5, (b) Run 18 showing the effect of setting /jj c = 0.95, (c) 
Run 19 showing the effect of using the LVK11 envelope binding energies, and (d) Run 35 showing the effect of increasing acE to 3.0. 
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Fig. 7. — Same as |Figurc 4| with each panel corresponding to the BSE runs shown in |Figurc 6| (a) Run 17, (b) Run 18, (c) Run 19, and 
(d) Run 35. 



and sdB + M dwarf binaries predicted by our BPS mod- 
els. The observed sample is highly incomplete so we can 
only test whether the BPS populations span the range of 
observed orbital periods and companion masses. Since 
we cannot compare the density of systems in the ob- 
served and synthesized sample, for each model we com- 
pute a model suitability parameter p, which is the frac- 
tion of observed systems that lie within boxes that en- 
close the modeled systems in companion mass-orbital pe- 
riod space. In some runs, the sdB + WD binaries exhibit 
a bimodal period distribution with a large gap between 
the short and long period systems. In cases where this 
gap in logP > 0.5 dex we draw separate boxes around 
the long and short period distributions. Furthermore, 
because the companion type is known for 26 binaries in 
the observed sample, we draw separate boxes around the 
sdB + WD and sdB + M dwarf binaries produced in our 
BPS runs and require that these observed systems fall 
within the appropriate box. The remaining 40 systems, 
for which the companion type is not known, are included 



in p if they lie in either box. We list p for each run in 
ITable II 

3.2. Parameter Variations 

To illustrate how parameter variations impact the pop- 
ulation of sdB binaries, we will move from one set of pa- 
rameters that uses the BSE defaults to another set that is 
comparable to that of run 2 in H03 with successive, cu- 
mulative changes. In all of these runs, we used the same 
common envelope ejection efficacy, acE = 0.75. Fur- 
thermore, the specific angular momentum carried away 
by mass lost from the system during stable mass trans- 
fer was also the same for each of these runs; the material 
lost during stable RLOF carried away the specific orbital 
angular momentum of the system (7 = 1). 

|Figurc 3 1 shows the population predicted by Run 1, 
which we will use as our fiducial run. The solid black 
lines bound the region populated by the H03 best fit 
models, and are shown for comparison. A single initial 
binary can produce multiple points on these diagrams 



Population Synthesis of sdBs 



11 



© 



© 




1 

log P [day] 



-1 1 

log P [day] 



Fig. 8. — Same as |Figurc 5| with each panel corresponding to the BSE runs shown in |Figure 6| (a) Run 17, (b) Run 18, (c) Run 19, and 
(d) Run 35. 



in proportion to the time it spends in the sdB "box" 
because of the resampling scheme used to generate a 
present day population. Sequences of symbols for the 
same initial system lie along that binary's evolutionary 
track and may be roughly horizontal if the companion 
evolves across the HG or vertical due to period evolution 
through magnetic braking or tidal exchange of spin and 
orbital angular momentum. We warn that the density 
of systems in a particular region cannot be read from 
these plots because several binaries are plotted on top 
of one anot her. Marginal distributions are discussed in 
|section 3.3| and the number of systems formed thorough 
each channel is given in |Tablc~T] 

In Run 1, almost all of the sdBs with hydrogen burning 
companions were formed through stable RLOF initiated 
while the primary was in the HG. The primaries in these 
system s had init ial masses M\q > 2 Mq, which can be 
seen in Figure 4 1 which shows the initial mass of the pri- 



mary Mio, initial period Pq, and mass ratio qo for the 
sdB 4- MS binaries. Their evolution was similar to that 
of Algol, the stable mass transfer eventually leading to 
mass ratio reversal and an increase in the orbital separa- 
tion. This resulted in sdBs with early type companions 
in systems with P > 10 days. In some cases, the com- 
panions have evolved into or across the Hertzsprung gap. 
In this run, there are nine sdBs that formed as the result 



of a CE phase that began while the primary was on the 
giant branch. These systems all have P < 1 day. 

When we allowed stars to ignite helium if their cores 
had reached a mass m c = 0.95m c (TRGB) before the en- 
velope was lost, far more sdBs were formed during CE 
evoluti on that be gan while the primary was on the RGB 
(Run 2, Figure 3| j) . The vast majority of the sdBs formed 
through CE evolution have G, K, or M type companions 
and the orbital periods of these binaries range from 0.05- 
220 days. |Figure 4p shows that primaries in these sys- 
tems (green) had initial masses in the range 0.8 — 2 Mq 
and that the initial periods of these systems were in the 
range 60-1100 days. In Run 1, these systems survived 
the CE phase, but the primary never ignited helium and 
became a helium WD without passing through the sdB 
"box." 

In Run 3, we have u sed the en velope binding energies 
computed by LVK11 (Figure 3;). These more tightly 
bound envelopes resulted in a much narrower range of 
orbital periods for post CE systems, none of which exceed 
10 days. Additionally, some systems were unable to eject 
the more tightly bound envelope and merged without 
ever producing an sdB. 

Settin g q CI jt — 1.5 on the RGB and 3.2 in the HG 
(Run 
CE phase 



Figure 3 1) allowed many systems to avoid the 
and undergo stable RLOF mass transfer on 
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Fig. 9. — Same as |Figurc 5| with each panel corresponding to the BSE runs shown in |Figure 6| (a) Run 17, (b) Run 18, (c) Run 19, and 
(d) Run 35. 



the RGB. These systems produce binaries consisting of 
a sdB with an A to G type MS companion and orbital 
periods in the range 55-660 days. As was the case with 
the systems formed by stable RLOF mass transfer ini- 
tiated in the HG, some of the companions have evolved 
into later type giants. There was also a slight decrease 
in the number of short period sdB + F, G, or K type 
dwarf systems because of the higher value of g C rit- With 
Qmt = 1-5, these systems underwent stable RLOF mass 
transfer that did not completely remove the envelope and 
never produced an sdB. Instead, they reach He ignition 
mass and detach from the Roche lobe before the entire 
envelope is transferred; they become horizontal branch 
stars. 

We considered a fixed fraction of mass loss during st a- 
ble RLOF mass transfer, /rlof = 0.5 (Run 5, Figure 3\i ) . 
This resulted in far more mass loss during stable RLOF 
mass transfer than in the previously discussed runs. This 
non-conservative mass transfer had two effects on the re- 
sulting population. First, the mass lost from the system 
carried away angular momentum which drove the bina- 
ries to shorter periods. Second, since less material was 
accreted by the secondary, the companions have lower 
masses and, thus, lower T e g than they did in the nearly 
conservative case. 

Finally, in |Figure 3f , we show the population predicted 



by Run 6, in which we used the default BSE values for 
the envelope binding energy. It is useful to compare t he 
post-CE populations of Run 6 and Run 2 (Figure 3d). 
Increasing q CI i t has depleted the high log P, high logT c ff 
corner of the post-CE population. 

Runs 1-6 produced a diverse set of populations of 
sdB + MS binaries. The orbital periods were in the 
range 0.04-663 days and companions were of every spec- 
tral type from B to M. Furthermore, it was possible 
to create these vastly differe nt populat ions while hold- 
ing ctcE and 7 constant. In Figure 5 we compare the 
sdB + WD and sdB + M dwarf binaries predicted by 
our BPS models with the observed sdB binaries from 



Morales-Rueda et al. (2003) and Copperwheat et al. 
( 201ip . 'lb facilitate the comparison with observed sdB 



binaries, we have computed a minimum mass for the 
WD or M dwarf companions: M 2 sin 3 i is "observed" for 
each theoretical present day system, via the mass func- 
tion f(M) = PKf/2irG = Ml sin 3 i/(Mi + M 2 ) 2 ; where 
K = sin i/P with i chosen so the distribution of 
cos i is uniform and the other quantities known from 
model output. Assuming M\ = 0.48 Mq then allows 
us to infer M 2 sin 3 i, and setting sini = 1 gives M m j n . 
Again, a single initial binary is responsible for several 
points on this diagram. In | Figure 5| the towers of points 
at nearly constant P correspond to a single, resampled 
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Fig. 10. — Same as |Figurc 3| for (a) Run 9 showing the effect of allowing the ejected material to carry away the specific angular momentum 
of the donor (7 = — 1), (b) Kun 10 showing the effect of setting /h c = 0.95, (c) Run 12 showing the effects of setting q cl -it = 1.5 and using 
the LVK11 envelope binding energies, and (d) Run 13 showing the effect of setting /rxof = 0-5. 



binary evolution track with an independent inclination 
for each point. In many of our BPS runs, one or two sdB 
+ neutron star binaries were formed. We included these 
binaries with the sdB + WD systems because the two 
are indistinguishable in terms of the simulated M m j n . 
The fraction of observed sy stems rep roduced by each 



Table 1 Clearly, Runs 1, 



run is listed in column (19) of [I 

3,4, and 5 are unable to reproduce many of the observed 
systems. The populations produced by Runs 2 and 6, on 
the other hand, are consistent with the observed popu- 
lation. However, despite the similarity amongst the sdB 
+ WD and sdB + M dwarf populations in Run 2, Run 
6, and the observed sample, the two models make strik- 
ingly different predictions for the popu latio n of sdB + 
MS binaries, as can be seen in |Figure~3p and |Figure 3f . 

In each run, 3 — 11% of the sdB + WD systems had 
previously been sdB + MS binaries. In these systems 
the primary becomes an sdB through an episode of sta- 
ble mass transfer. Eventually this sdB evolves into a WD 
and after the secondary begins to ascend the RGB the 
system undergoes CE evolution resulting in a sdB + WD 
system. The period and mass distributions of the sdB 
+ WD binaries that form through this channel overlap 



with those of the sdB + WD binaries formed in systems 
that have not previously undergone an sdB phase, which 
makes it difficult to distinguish between the two subpop- 
ulations. 

3.2.1. Influence of acE 

We demonstrate the effect of varying acE m Figures 
5-8, again with succ essive, cumulative changes to the pa- 
rameters. Figure 6 1 shows Run 17, which was the same 
as our fiducial run, shown in |Figure 3 1, except we in- 



creased acE to 1.5. In this case, the orbital energy was 
efficiently transferred to the CE as the stars spiraled to- 
wards one another and many more systems were able to 
eject the CE before merging. The result is a band of 
systems with P < 5 days with companions spanning all 
types from early A to M. |Figure 7 % shows that these 
post-CE systems occupy a different region of parameter 
space than most of the post-CE systems seen in Runs 
1-6, namely that M±o > 2 Mq for these systems. The 
primaries in these binaries were able to ignite helium 
non-degenerately in their cores, despite the loss of the H 
envelope in the CE phase. 
For Run 18, we kept acE = 1-5 and set /h c = 0.95. 
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Fig. 11. — Same as |Figure 4| with each panel corresponding to the BSE runs shown in |Figure 10] (a) Run 9, (b) Run 10, (c) Run 12, and 
(d) Run 13. 



This resulted in a large number of additional post-CE 
systems, most of which had longer periods than the post- 
CE systems produced in Run 17. These systems have 
longer periods because their formation required that the 
primary did not fill its Roche lobe until it was near the 
TRGB, otherwise the CE phase would have begun prior 
to the core reaching the mass required for helium igni- 
tion. This, in turn, required binaries with large initial 
separations, in order to incre ase the siz e of the Roche 
lobe. Co mparing pa nel (a) of Figure 4 with panels (a) 
and (b) of Figure 7 one can clearly see the different re- 
gions of initial parameter space that produced post-CE 
systems as a result of increasing q;ce (the band of sys- 
tems with Mio > 2 Mq) and as a result of reducing fa e 
(the systems with M\o < 2 Mq). 

We used t he LVK11 binding energies in Run 19, shown 
in |Figure 6} :. Even though the orbital energy was effi- 
ciently transferred to the envelope in this run, the higher 
envelope binding energies in Run 19 resulted in most 
of the systems being driven to merger, which in Run 



18 ignited helium non-degenerately and survived the CE 
phase. In [Figure Oh 1 we show Run 35, which is identical 
to Run 19 except we have set q;ce = 3.0. Allowing the or- 
bital energy to be transferred to the envelope even more 
efficiently allows systems to avoid mergers, even with the 
i ncreased b inding energies of LVK11. 

[Figure 9" shows that Run 35 is the only one of these 
four models that is consistent with observations. The 
short period binaries in Run 18 do span the observed 
range of period and mass, however the large gap in the 
sdB + WD period distribution is problematic. Several 
of the observed sdB binaries known to have WD com- 
panions fall within the gap, which makes the population 
inconsistent with the sample. 

3.2.2. Influence of 7 

Figures 9-12 illustrate how changing the specific an- 
gular momentum carried away by material lost from the 
system during mass transfer affects the population of sdB 
+ MS binaries. In each of the runs shown in these fig- 
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Fig. 12. — Same aslFigurc 5|with each panel corresponding to the BSE runs shown inlFigurc 10| (a) Run 9, (b) Run 10, (c) Run 12, and 
(d) Run 13. 



ures, the material lost from the system carries away the 
specific angular momentum of the donor et^fiorb- The 
angular momentum lost per unit mass in these runs is at 
least a factor of four lower (the primary is the donor so 
a jay = (g+1) an d q > 1) than in the previously discussed 
runs. |Figure lOfc , shows Run 9, which used the same pa- 
rameters as our fiducial run except for angular momen- 
tum loss. Reducing the amount of angular momentum 
carried away by ejected material resulted in the forma- 
tion of sdB + A or F dwarf binaries with P > 10 days. 
These binaries were not produced in Run 1 because an- 
gular momentum loss during stable mass transfer drove 
these systems into contact, resulting in a merger. Ad- 
ditionally, there are some sdB + A or F dwarf binaries 
formed after a CE phase that begins while the primary is 
in the HG. These systems began CE evolution while the 
stars were more widely separated than the same systems 
in Run 1, a result of the diminished angular momentum 
lo ss during R L O F . 

|Figurc Kip shows the population produced by Run 10, 
where we have allowed stars to ignite helium shy of the 
TRGB. The post-CE population produced as a result is 
nearly identical to that produced in Run 2 because these 
systems did not undergo stable mass transfer and are 
therefore unaffected by the altered angular momentum 
loss condition. We have plotted Run 12 in |Figure lOp . 



In this run we have changed both <7 cr it and t he the bin d- 
ing energy. By comparing this run to Run 4 (Figure 31), 



one can see that reducing the angular momentum carried 
away by ejected material has resulted in a slight decrease 
in the number of sdB + G dwarf systems with P > 10 
days. Without the angular momentum loss, mass trans- 
fer ended in these systems before the primary's envelope 
was completely stripped. Finally, in Figure 10 i, we show 
Run 13 and the effect of forcing non-conservative mass 
transf er by setting /rlof = 0.5. As with Run 6 (|Fig- 



ure 3 ? ) , T e g is lower for the companions in this case be 



cause they have not accreted as much material and are 
therefore less massive. However, since the lost material is 
removing far less angular momentum from the system in 
Run 13 than in Run 6, the orbital periods of the binaries 
produced by stable RLOF mass transfer in Run 13 are 
only slightly reduced, with most remaining longer than 
10 days. 

Comparing Figure 11 (7 = —1) with Figure 4 (7 = 1) 
reveals that the amount of specific angular momentum 
carried away by ejected material only weakly impacts 
which initial binaries will become sdB binaries. While 
these runs were useful to illustrate how changing 7 affects 
the population of sdB + MS binaries produced in our 
BPS models, Figure 12 shows that many of these runs 



do not reproduce the observed systems. 
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Fig. 13. — Histograms of the companion logT c fj distribution 
for the runs shown in Figures 1-11. Both true sdB and post-RGB 
binaries are included. 

3.3. Distributions 

Another useful way to investigate the BPS populations 
is to lo ok at the distributions of their properties. |Fig-| 
|ure 13| shows the distribution of the companion's log T c ff 
in each of the runs discussed above. These distributions 
are not readily apparent in Figures 3, 6, and 9 due to 
severe overplotting in some regions. The bin-size is 0.1 
in log T e ff, with bin edges centered at multiples of 0.1 
in log T c g. There are two distinct groups, one at high 
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Fig. 14. — Histograms of the age of the sdB + MS and post-RGB 
-I- MS binaries for the runs shown in Figures 1-11. 



log T e g corresponding to the HG RLOF systems and one 
at low log T e g due to the RGB CE II systems. The for- 
mer group is always present. The RGB CE II group is 
dominant when present, but only appears in runs with 
/hc = 0.95. The region between these groups can be 
populated by systems formed through either the RGB 
CE I channel or the RGB RLOF channel. When the 
RGB CE II group is not present, the post-RGB systems 
dominate the distribution for T c g < 3.85. The same 
grouping, by formation channel, is seen in the age dis- 
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Fig. 15. — Histograms of Af s( jB f° r the runs shown in Figures 
1-11. Both true sdB and post-RGB binaries are included. 



tributions of these runs, which are shown in Figure 14) 
Most of the post-RGB binaries are formed in binaries 
with Mio < 2M and they produce the distinct feature 
at high age seen in Runs 1, 17, and 9. 

The m asses of the sdBs in the modeled populations are 
shown in Figure 15 When /jj e = 0.95, the distribution 
peaks near the canonical value of 0.48 M . The majority 
of sdBs in these runs are formed through the RGB CE II 
channel by progenitors that ignite He degenerately near 
the TRGB. Stars that ignite He non-degenerately on the 
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Fig. 16. — Histograms of the orbital periods of helium burning 
sdB + early F binaries (7350 > T cft > 6700, hatched) from the 
"focused" BPS runs and post-RGB + early F binaries (red) from 
the twelve runs with a model suitability parameter p > 0.9 (see 
|scction 3.1| . The post-RGB histograms have been rescaled to be 
m proper proportion to the "focused" models and we have fixed 
the larger of the two density distributions maxima to unity, to em- 
phasize the relative differences between distributions on the same 
scale. 
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RGB produce the sdBs that populate both the low and 
high mass tails of these distributions. Primaries with 
Mio ~ 2 Mq that begin transferring mass in the HG can 
become low mass sdBs with 0.32 < M s( jb ;$ 0.4. Interest- 
ingly, our models still produce high mass sdBs, despite 
the fact that we do not follow mergers. More massive 
(A^io ^ 4 Mq) primaries undergoing stable RLOF mass 
transfer that begins in the HG yield sdBs with masses 
as high as 0.79 Mq. In each run, the post-RGB stars 
populate a separate, low-mass region. 

4. CONSTRAINTS ON BPS PARAMETERS 



In |section 3.1| we defined the run suitability parame- 
ter for each run, p, that is the fraction of observed short 
period sdB binaries spanned by the short period sdB + 
WD and/or sdB + M dwarf binaries in each BPS model. 
Twelve of the parameters sets explored here were able 
to reproduce more than 90% of systems in the observed 
sample, despite making very different predictions about 
the sdB + MS population. For any value of any one 
parameter, a combination of the others exists which is 
qualitatively consistent with the observed sample. In 
fact, the only parameter varied in our study that was well 
constrained by current observations is /h c - Only mod- 
els with /ho = 0.95 produced sdB + WD and sdB + M 
dwarf systems with P > 1 day, which are required by ob- 
servations. However, this is not necessarily a statement 
about the minimum core mass at which helium ignition 
occurs. With /h c = 1-0 a primary that ignites helium de- 
generately (i.e., one with Mi < 2 Mq) could only form 
an sdB if it filled its Roche lobe at exactly the same time 
it reached the TRGB and ignited helium. This scenario 
depends on many assumptions built into BSE about sin- 
gle star evolution, the size of the Roche lobe, and the 
onset and duration of unstable mass transfer, such that 
there is a narrow range of qo and Pq with a given Mio for 
which it is possible, as was pointed out by |Mengel et al.| 
(1976). (Compare panels a and b in Figures 4 and 5 for 
example). In practical terms, changing /h extends the 
range of initial orbital periods for which a system with 
a primary of M < 2 Mq and q > q CY i t will produce a 
binary containing an sdB, thus adding some flexibility to 
the BSE algorithms. This flexibility is essential to repro- 
ducing the observed population, but it is impossible to 
determine exactly which physical process or combination 
of processes requires this freedom. Further observations 
are needed to better constrain binary evolution parame- 
ters and go beyond the limits of existing studies. 

Still, twelve combinations of the remaining parame- 
ters produced a population consistent with the observed 
sample. Models that used the LVK11 binding energies 
required values of acE > 0.75 to produce binaries similar 
the observed systems with the longest periods and high- 
est companion masses. Nevertheless, fixing the value of 
any one parameter, we can find a set of values for the 
remaining parameters that produces a population con- 
sistent with the observed sample. Note, however, that 
if the BPS parameters are correlated, this may not be 
possible. We conclude that the short period sdB + WD 
and sdB + M dwarf systems that constitute the presently 
well-explored sample are fairly insensitive to the param- 
eterizations used in our BPS models, and propose that 
observations of sdB + early F dwarf binaries could better 
them. 



4.1. The Advantage of sdB + Early F dwarf Binaries 

Our models suggest that the sdB + early F dwarf 
(7350 > T ff > 6700) binary population can discrimi- 
nate between different binary evolution scenarios. Stars 
in this temperature range are cool enough that the sdB 
will produce a measurable far-ultraviolet excess, without 
being so cool as to make selection based on optical/IR 
colors difficult. We show normalized histograms of the 
helium burning sdB + early F binary period dist ributions 
for the twelve runs with p > 0.9 in Figure 16 (hatched 
histograms). To ensure that our models properly sam- 
pled the period distribution of sdB + early F dwarf bina- 
ries, we ran additional, "focused" BPS models for each 
of these twelve runs. We drew 10 6 initial binaries from 
the same distributions described above, but we restricted 
the range of Mio, a, and qo to maximize the number of 
sdB + early F dwarf binaries. These runs produced 200- 
1100 ti mes as ma ny sdB + early F dwarf binaries as those 
listed in |Tablc Tj We find good agreement between period 
distributions in these models and those described above 
and we plot the distributions for the "focused" models. 

These twelve sets of parameters predict populations of 
sdB + early F dwarf binaries that are useful in constrain- 
ing <7crit> 7, and acE- First, these systems will only have 
P > 75 days if they are formed through stable RLOF 
mass transfer that begins while the primary is on the 
RGB, sec Figure 16 panels b, c, d, g, h, k, and 1. This 
scenario is only possible if q cl -i t > 1, which is almost 
never the case given Equation 2| The presence of sdB 
+ F dwarf binaries with P > 75 days would be evidence 
that mass transfer beginning on the RGB can be stable. 
Furthermore, if the mass ejected from the binary dur- 
ing non-conservative mass transfer only carries away the 
specific angular momentum of the primary (7 = —1), the 
orbital periods of the sdB + early F bi naries would be 
even longer, P > 250 day, see |Figure 16| panels d, k, and 
1. 

On the other hand, the absence of long period (P > 75 
day) sdB + F dwarf binaries would indicate that mass 
transfer beginning on the RGB is unstable and that such 
systems cannot avoid CE evolution. If systems that be- 
gin mass transfer while the primary is on the RGB do 
undergo CE evolution, they will only survive as short pe- 
rio d (P < 5 d ays) sdB + F d warf binari es if «ce> 0.75, 
Figure 16 panels c-1. From [Figure 16 a, it appears that 



there are short period sdB + F dwarf binaries that sur- 
vive with acE = 0.75, however the figure is misleading 
because there is only a s ingle sdB + early F dwarf binary 
produced in Run 2 (see|Figure 3 3). 

In panel j of |Figurc 16| we see a population of sdB + F 
binaries with P > 10 days. These binaries form through 
stable RLOF mass transfer that begins while the primary 
is in the HG, but only if material ejected from the sys- 
tem carries away the specific a ngular mom entum of the 
primary (see panels a and b of Figure 10). These bina- 
ries should be distinguishable from those formed through 
stable mass transfer beginning on the RGB because they 
all have P < 70 day. The constraints placed on binary 
evolution scenarios by each possible population of helium 
burning sdB + F early dwarf binaries can be summarized 
as follows: 

1. No sdB + early F binaries: unstable mass transfer 
on the RGB and acE ^$ 0.75 
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2. Only sdB + early F binaries with P > 75 days: 
stable mass transfer on the RGB, «ce ^ 0.75 

a) If sdB + early F binaries have P > 250 day 
then 7 = — 1, otherwise 7 = 1 

3. Only sdB + early F binaries with P < 75: unstable 
mass transfer on RGB and acE > 0.75 

a) If sdB + early F binaries have P > 10 day 
then 7 = — 1, otherwise 7 = 1 

4. Long and short period sdB + F: stable mass trans- 
fer on the RGB and acE J$ 0.75 

a) If sdB + early F binaries have P > 250 day 
then 7 = — 1, otherwise 7 = 1 

Clearly, the helium burning sdB + early F dwarf binary 
population has exceptional diagnostic power. However, 
such binaries are largely absent from catalogs because 
sdBs are discovered in optical color surveys that search 
for faint blue stars; sdB + early F binaries may be 
misse d because the lu minous companion outshines the 
sdB ( |Wade et al.pOlOl ). 

There is an additional complication to studying the 
sdB + F binaries. Many of the post-RGB objects have F 
type companions (see Figures 3, 6, and 9), and we plot 
their period distribution (rescaled to be consistent with 
the "focused" BPS models) as the red histogram in Fig- 
This overlapping population requires that the 
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nature of the F star's companion be carefully confirmed 
before the system can be used to constrain binary evolu- 
tion models, because the period distribution of post-RGB 
binaries is not as sensitive to the parameters explored 
here. 

5. COMPARISON WITH PREVIOUS STUDIES 

In Figures 3, 6, and 9 we have indicated the region in 
log P — log T e g space occupied by the binaries predicted 
in the best fit model chosen in H03 (see their Figure 15). 
In our Run 6, we used parameters similar to H03 and 
were able to reproduce many features of the binary pop- 
ulation presented there. The orbital periods predicted in 
our model range from 1 hr to 663 days which is similar to 
the 0.5 hr to 500 d ay range pr e dicted in H03. Although 



the calculations of |Han e t al 
be possible to create sdB 



(2002) show that it may 
MSbmaries with P > 1000 



days, neither our BPS models nor the H03 BPS models 
produce such systems. One notable difference is that our 
model predicted systematically lower temperatures for 
the companions in systems that formed through stable 
RLOF mass transfer. This difference can be accounted 
for, in part, by the different approaches used in the BPS 
calculations. H03 assumed that RLOF mass transfer was 
instantaneous, and appli ed the resul t s of d etailed binary 
calculations presented in Han et al. (2002) to determine 
the minimum core mass required at the onset of RLOF 
to produce an sdB. If RLOF began when the primary's 
core was massive enough (O.28A/ < m c < 0.39Af Q , de- 
pending on the star's total mass), the entire envelope was 
stripped and half of the envelope mass was transferred 
to the companion (Z. Han, private communication). Be- 
cause the core was not allowed to grow during RLOF, this 
resulted in less massive sdBs and excess material being 



transferred to the companions. Our BSE models allowed 
the mass donor to evolve during RLOF which decreases 
the mass of the envelope and reduces the amount of ma- 
terial transferred to the companion by ~ 0.1 M . The 
lower temperatures predicted in our models result in a 
smaller "gap" in the temperature distribution for sys- 
tems with P < 10 days. In H03 there are no systems 
with 3.8 < logTeff < 4.0 at periods shorter than 10 days. 
In Run 6, the gap only extends from log T e g = 3.8 — 3.85. 

The nature of the post-CE systems in some of our mod- 
els differed g reatly from the predictions of H03 an d Nclc 



mans 



(2010). In the models that take g cr it from Equa- 



tion 2 the periods of the post-CE binaries extend to 



much higher values than those predicted by H03. The 
maximum period of these post-CE systems ranges from 
181 days with acE = 0.75 to 543 days with acE = 3.0 
(e.g., Run 34). These long period post-CE systems are 
produced by binaries with wide initial separations and G 
star p rimaries. The BPS models presented in |Nelemans| 
(20101 use a different prescription for the CE ejection 



process in which it is possible for the period of a system 
to increase during during the CE phase. This results in 
post-CE binaries with P > 1000 days, far longer than 
the long period systems predicted in our models. The 
short period, post-CE systems with G or later compan- 
ions predicted by our models are consistent with those 
predicted by H03. This is contradicted by the lack of 
observed short period sdB + G dwarf systems. However, 
some of our models show there are far more sdB + G 
dwarf systems with long periods than with short, w hich 
is consistent with the findings of |Green et al. (2001 1 and 
Copperwheat et al. (20111. 



6. CONCLUSION 

We have presented a set of BPS models that explored 
the population of sdBs in binaries with hydrogen burn- 
ing companions. We found that this population is highly 
sensitive to the parameters used to simulate binary evo- 
lution. Plausible variations in the assumptions about 
the minimum core mass for helium ignition, the enve- 
lope binding energy, the common envelope ejection effi- 
ciency, the amount of mass and angular momentum lost 
during stable mass transfer, and the criteria for stable 
mass transfer result in populations of sdB + MS binaries 
with drastically different orbital period and companion 
temperature distributions. Our BSE models also suggest 
that the population of sdB + WD and sdB + M dwarf 
binaries used to constrain the values of these parameters 
cannot do so unambiguously. Twelve very distinct pa- 
rameter sets can reproduce over 90% of the systems in 
the currently well-studied sample. 

The period distribution of helium burning sdB + early 
F dwarf binaries, on the other hand, can differentiate be- 
tween the binary evolution scenarios described by these 
parameterizations. The presence of sdB + early F dwarf 
binaries with P < 5 days would indicate that orbital en- 
ergy is efficiently transferred to the CE. Additionally, the 
existence of sdB + early F dwarf binaries with P > 75 
day would require that mass transfer beginning while 
the primary is on the RGB is stable. Details of the sdB 
+ early F dwarf binary distribution also probe angular 
momentum loss during non-conservative mass transfer. 
Additionally, the models predict a population of binaries 
containing post-RGB objects that lie in the sdB "box" 
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that is, like the population of sdB + early F dwarf bina- 
ries, largely absent from catalogs. 
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0.00 





00 


91.23 


0.48 


96.66 


1.41 


1.93 





88 


m 

a 


28 


1 


50 


-1 


BSE 


1 . 5 


■3.2 


LVK10 


0.95 


11401 


ZU4ool 


3.87 


5.53 





(JO 


84.28 


0.46 


94.13 


4.83 


1.03 





78 


CO 


29 


1 


50 


-1 


0.5 


1.5 


3.2 


LVK10 


0.95 


12115 


f)AO A 77 

208477 


4.18 


5.70 





00 


82.62 


0.85 


93.35 


4.24 


2.42 





85 


cr. 


30 


1 


50 


-1 


BSE 


1.5 


3.2 


BSE 


0.95 


13220 


22o22o 


3.82 


5.21 





00 


82.56 


0.78 


92.38 


5.43 


2.20 





91 




31 


1 


50 


-1 


0.5 


1 . 5 


■3.2 


BSE 


1.00 


4207 


34600 


25.19 


14.12 





00 


24.55 


5.17 


69.03 


16.49 


14.48 





26 




32 


1 


50 


-1 


0.5 


1 . 5 


■3.2 


LVK10 


1.00 


3204 


22803 


38.22 


21.43 





00 


0.03 


7.74 


67.42 


10.50 


22.08 





08 




33 


3 


00 


1 


BSE 


BSE 


4.0 


BSE 


1.00 


8169 


o o r- C O 1 

Z3o5z4 


1.34 


0.00 





00 


92.31 


0.39 


94.35 


5.46 


0.19 





71 




34 


3 


00 


1 


BSE 


BSE 


4.0 


BSE 


0.95 


16979 


430671 


0.73 


0.00 





00 


95.39 


0.22 


96.52 


3.38 


0.10 





98 




35 


3 


00 


1 


BSE 


BSE 


4.0 


LVK10 


0.95 


11161 


234647 


1.33 


0.00 





00 


95.85 


0.37 


97.56 


2.23 


0.21 





92 




36 


3 


00 


1 


BSE 


1.5 


3.2 


LVK10 


0.95 


12128 


O/ ■ C\C\ O A 

260984 


1.54 


4.44 





00 


82.43 


0.33 


88.74 


11.01 


0.25 





89 




37 


3 


00 


1 


0.5 


1.5 


3.2 


LVK10 


0.95 


11913 


246390 


1.58 


4.63 





00 


87.31 


0.07 


93.59 


5.79 


0.62 





91 




38 


3 


00 


1 


0.5 


1 . 5 


3.2 


BSE 


0.95 


17611 


433900 


0.91 


2.63 


2 


63 


92.00 


0.05 


95.65 


4.00 


0.34 





98 




39 


3 


00 


1 


0.5 


1.5 


3.2 


BSE 


1.00 


8753 


239108 


1.65 


1.89 


1 


89 


89.75 


0.10 


93.50 


5.88 


0.62 





75 




10 


3 


00 


1 


0.5 


1.5 


3.2 


LVK10 


1.00 


3083 


53220 


7.33 


8.51 





00 


59.40 


0.33 


75.57 


21.57 


2.87 





16 




41 


3 


00 


-1 


BSE 


BSE 


4.0 


BSE 


1.00 


10007 


251117 


3.89 


0.00 





00 


86.94 


0.42 


92.44 


6.08 


1.48 





77 




42 


3 


00 


-1 


BSE 


BSE 


4.0 


BSE 


0.95 


18818 


445368 


2.19 


0.00 





00 


92.24 


0.23 


95.34 


3.82 


0.84 


1 


00 




43 


3 


00 


-1 


BSE 


BSE 


4.0 


LVK10 


0.95 


12932 


245133 


3.98 


0.00 





00 


91.75 


0.40 


96.14 


2.33 


1.54 





92 




44 


3 


00 


-1 


BSE 


1.5 


3.2 


LVK10 


0.95 


13066 


267069 


2.96 


4.23 





00 


80.55 


0.37 


88.11 


11.11 


0.78 





88 




45 


3 


00 


-1 


0.5 


1.5 


3.2 


LVK10 


0.95 


13607 


256735 


3.39 


4.63 





00 


83.79 


0.70 


92.52 


5.53 


1.95 





91 




46 


3 


00 


-1 


BSE 


1.5 


3.2 


BSE 


0.95 


19373 


445589 


1.96 


2.67 


2 


67 


89.59 


0.42 


94.73 


4.15 


1.11 





98 




47 


3 


00 


-1 


0.5 


1.5 


3.2 


BSE 


1.00 


10321 


246973 


3.53 


1.98 


1 


98 


86.89 


0.75 


93.33 


4.66 


2.01 





72 




18 


3 


00 


-1 


0.5 


1.5 


3.2 


LVK10 


1.00 


4581 


59674 


14.60 


8.19 





00 


52.97 


3.01 


78.78 


12.82 


8.39 





43 





a Columns (11) - (18) are percentages of the number of systems given in column (10) formed by each channel. The values in columns (11)-(15) sum to the 
value in column (1 6). 



'See 
c See 



Equation 3 



Equation z 



d This run most closely resembles the best fit run in H03. 



